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O ' Abstract 



We study the time evolution of a qubit linearly coupled with a quantum environment under a sequence 
of short 7r pulses. Our attention is focused on the case where qubit-environment interactions induce 
the decoherence with population decay. We assume that the environment consists of a set of bosonic 
excitations. The time evolution of the reduced density matrix for the qubit is calculated in the presence 
of periodic short tt pulses. We confirm that the decoherence is suppressed if the pulse interval is shorter 
than the correlation time for qubit-environment interactions. 
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1 Introduction 

To employ a two-level quantum system as a qubit for quantum computations, phase coherence must be well 
preserved in it for a long time. However, since any quantum systems interact with their environment, deco- 
herence inevitably arises more or less. Thus, several schemes to suppress decoherence have been proposed. 1-4 
Among them, the decoherence suppression scheme which uses a sequence of short tt pulses 3-7 is a promising 
O-i, candidate for actual applications since it is relatively easy to implement in experimental situations and can 
be readily applied to multi-qubit systems. We call this scheme the pulse control of decoherence. The pulse 
control has been proposed by Ban, 3 and by Viola and Lloyd. 4 They considered a qubit under a sequence 
of short tt pulses assuming that qubit-environment interactions induce the decoherence without population 
decay (i.e., dephasing). They calculated the time evolution of the qubit, and showed that the dephasing 
can be suppressed by applying periodic tt pulses if the pulse interval is shorter than the correlation time for 
^ ■ qubit-environment interactions. It has been proposed that we can even suppress the decoherence with pop- 
ulation decay by applying tt pulses. 5 Let us consider a qubit linearly interacting with a boson environment. 
The Hamiltonian is Tto = Tts + ~Hb + H-sb with 

H s = ^oy^lB, (1) 
H B = 1 S <8>X^ 6 K- ( 2 ) 

k 

Here, Ti.s and TL-q describe the qubit and the environment, respectively, and 7^sb describes their mutual 
interactions. We have used the pseudo-spin representation in expressing the qubit. Ban, and Viola and 
Lloyd, considered the case where TLsb is given by 

W S b = ^h(gla z ®b\ + g k a z ®6 fc J , (3) 



which induces the dephasing. They proposed that the dephasing can be suppressed by applying periodic 
tt pulses around the x-axis in the pseudo-spin space. Let p(t) be the density matrix for the whole system 
consisting of the qubit and the environment. We can simply integrate the Liouville-von Neumann equation 
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for p(t) under the ir pulses since [7^s, Wsb] = 0, and explicitly examine the effectiveness of the pulse control. 
In contrast, if 

n SB = ^2h (g* k <T- ®b\ + g k a+ <g> 6*) , (4) 

k 

the decoherence with population decay arises. It has been proposed that the decoherence with population 
decay can be suppressed by a sequence of ir pulses around the z-axis. 5 However, its effectiveness has not 
been quantitatively examined so far. 

Generally, a qubit experiences not only the dephasing but also the decoherence with population decay. 
Thus, it is desirable to examine the effectiveness of the pulse control in more general situations beyond the 
case treated in refs. 3 and 4. In this paper, we present a formulation to describe the time evolution of the 
qubit coupled with a quantum environment under a sequence of tt pulses. Our formula is based on the time- 
convolutionless (TCL) projection operator approach, and is applicable to the case where the decoherence 
with population decay is present. An extension to more general cases where both the dephasing and the 
decoherence with population decay play a role is straightforward. 



2 Formulation 

Let us consider the case where the qubit-environment interaction is given by eq. (jlj and periodic tt pulses 
with an interval At are applied. Let 5t be the duration of each pulse. The influence of the tt pulses is 
described by H p (t) = V(t)a z (g> 1b, where 

r for m(At + St) < t < m(At + St) + At, 

V(t) = ( irh for + + At < t < (m + 1)1 At + St) ® 

v 2ot 

with m = 0, 1, 2, • • • . Thus, the total Hamiltonian is TL{t) = TCq + TC p (t). We decompose TL(t) as H(t) = 
TLi{t) + TtsBi where TLi{t) = 7is + ~Hb + TC p (t). The whole system is described by the density matrix p(t), 
which obeys the Liouville-von Neumann equation 

J^p(i) = L(t)p(t) (6) 



with 

L(t) ■ = - 

h 



L{t) . =-I[W(t),-]. (7) 



It is convenient to introduce p(t) which is defined as 

p(t)=%e-fi* rLl Wp(t), (8) 
where T & is the anti-time-ordering operator and L\{t)- = — j[[Hi(t) , •]■ From eq. ©, we obtain 

^- t p(t) = L S B(t)p(t), (9) 
where L S B(t)- = -^[H-sb^), ■}. Here, 7Y S b(*) is given by 

n S B(t) = T a e^ /o^it^WsBTe-s/o^iW, (10) 

where T is the time-ordering operator. We derive an explicit expression of TCsB^t)- Noting that 7i\{t) = 
H s +Hb+ H p (t) and [Hs + H B ,H p (t)] = 0, we obtain 

n SB (t) = T^^ dTH ^HsB{t)Te-TJo dTH ^), (11) 

where 

n SB (t) = es( w s+w B )t WgBe -^(H s +H B )t 

= J2 H ® 6^*""°)* + g k cy + ® b k e-^ k -^ . (12) 
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We here take the limit of 5t — > 0. This approximation is justified as long as 5t is much shorter than the 
typical time scales for the decoherence. In the interval between the iVpth pulse and the (iV p + l)th pulse 
(i.e., N p At <t<(N p + l)At ), we find that 



i /"* 7T 

- / cItH p (t) = i-N p a z 1 B . 
"jo ^ 



(13) 



We thus obtain 



which results in 



= (-l) JV ^ SB (t), (14) 

Lss(t)- =-(-l) N ^[HsB(t), •]. (15) 

Note that each tt pulse changes the sign of Lsb(£). Due to the periodic sign changes, the qubit-environment 
interaction is effectively reduced. This is the reason why the decoherence is suppressed by the tt pulses. 

Our interest is focused on the time evolution of the qubit, which is described by the reduced density 
matrix defined by 

(16) 

where pn (poo) represents the population in the upper (lower) state and the off-diagonal terms characterize 
the coherence of the qubit. In deriving the equation of motion for ps(t), we assume that the qubit and the 
environment are uncorrelated at t = 0, i.e., 

p(0) =Ps(O)0p B (O), (17) 
and the environment is initially in thermal equilibrium at temperature T, i.e., 

p B (0) = p B = IJ f 1 - e"^ fc ) e~^ kb i bk (18) 

k 

with j3 = 1/ksT (ks- the Boltzmann constant). We introduce the projection operator P which is defined as 
P- = TrB{-} <8> pb- We observe that P p(t) = ps(t) S3 pb- Adapting the TCL formalism 8 with P to eq. ©, 
we can derive the equation of motion for ps(t), 

-p s (t) = Tr B ^K 1 (t)Pp(t), (19) 

3=1 



t rt 



•i r l j-2 



where 

Kj(t) = I dh I dt 2 --- I dt j -i{LsB(t)LsB{ti)---L S B(tj-i)) .c- (20) 
Jo Jo Jo 

Here, (• • • ) ox . denotes an ordered cumulant 8 
(L SB (t) ■ ■ ■ £sb(*;-i))o.c. = ^(-l)^ 1 II pL SB(t) ■ ■ ■ L SB (t ql )PL SB (t q 2) ■ ■ ■ L S B(t q3 )PL SB (t q4 ) ■■■ . (21) 

The sum in eq. 1)21(1 should be taken over all possible divisions by P keeping the chronological order specified 
by t > ■ ■ ■ > t q i, t q 2 > ■ ■ ■ > t g 3, ■ ■ ■ and so on, and M is the number of P in each term. 

Note that Kj(t)Pp(t) with an odd j vanishes in our case. The qubit-environment interaction is not 
strong in actual situations, so we are allowed to retain only a few lower-order terms in the right-hand side 
of eq. I|19|) . For simplicity, we here approximately neglect the higher-order terms Kj(t) with j > 4. The 
role of such higher-order terms will be discussed elsewhere. We thus obtain 

^ps(i) = Tr B K 2 (t)Pp(t) 

= Tr B {Z SB (i)Wti){ps(t)®PB}}. (22) 
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Using eqs. (JUJ) and (15)1. we obtain for N p At < t < {N p + l)Ai, 





dt hi{t) = -ln(t)pu(t) + Vii(t), (23) 
§- t Pw(t) = -7io(t)pio(0, (24) 



where 



7ii (t) = Vbfcl 2 (2n B (w fe ) + l) / di 1 2cos(o; fe -u;o)(t-ti), (25) 

7lo(t) = V| 5fc | 2 (2n B (^) + l) / rftie 1 ^-^^*-*!), ( 26 ) 

lu(i) = Y^bfc| 2 "-B(^fc) / dii 2cos(u fc - uj ){t - h), (27) 
k J "p 

r ft r(N P -j)At 

dt x --- = / dh V (-iy / dh (28) 

JN P J NpAt j=Q J(N p -l-j)At 

and n B {u) is the Bose-Einstein distribution function. In deriving eq. (|23j) . we have used the relation of 
Pn(t) + Poo(t) = 1. Introducing the spectral function I{uj) defined by 

= Y,\9k\ 2 5(LU-LU k ), (29) 

k 

we rewrite eq. (|25j) as 

ln (t) = J dujI{uj)[2n B (uj) + l) J dh2 cos{uj - LO ){t -h) (30) 
and eqs. (|26|) and (|27j) are rewritten in similar manner. 

3 Results 

On the basis of the resulting equations, we study the time evolution of the reduced density matrix ps(t) in 

the presence of periodic tt pulses with an interval At. We numerically calculate p~u(t) and pw{t) assuming 
that the initial state is (|0) + |l))/\/2, i.e., 

*<"> - m) ■ < 31 > 

and that the spectral function is given by 

I(u)=ue- U ^, (32) 

where u c is cut-off frequency. The following parameters are employed: u> c /^o = 5.0 and feT/(fco) = 0.1. 

Then, the correlation time t c = 2ir/u; c for qubit-environment interactions is obtained as i c /(27r/u;o) = 0.2. 
If the 7T pulses are absent, pn(t) approaches to 

e -rkuo/(2k B T) 

^°°) = e+ ^o/(2k B T) + e -Hu / { 2k B T) ~ °( 10_5 ) ( 33 ) 

with increasing t. We treat the cases of At / (2tt /uq) = 0.016 and 0.032, for which the condition At <C i c 
holds. We display pu(t) and |/5io(i)| in Fig. ^a) and^b), respectively. Note that the decrease of p~\\ (|/5io|) 
represents the population decay (decoherence). From Fig. ^a), we observe that the population decay is 
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Figure 1: Time evolution of (a) p\\ and (b) |pio| at U-qT = 0.1 x hjjQ when the initial state is (|0) + 
|l))/v2- The dashed lines and dash-dotted lines correspond to the cases of /S.t / (2it / ujq) = 0.016 and 0.032, 
respectively. The solid lines represent the case where ir pulses are absent. 



suppressed by the ir pulses and the suppression becomes notable with decreasing At. Figure ^b) shows 
that the coherence of the qubit is recovered by the tt pulses. 

We apply the tt pulses around the z-axis, so each pulse does not directly modify the population. Note 
that reduction of the population decay is attributed to the periodic sign changes of the qubit-environment 
interaction caused by the tt pulses. Here we mention to an another way to reduce the population decay. 
The population decay can be suppressed without using tt pulses if we simply raise T. Indeed, since we have 
set pn(0) = poo(0) = 0.5, the population decay disappears in the extreme limit of k&T 3> friOQ. However, 
this inevitably accelerates the decoherence. 

4 Summary 

We have presented a formulation to describe the time evolution of a qubit coupled with a quantum envi- 
ronment under a sequence of short tt pulses. We have applied it to the case where the environment induces 
the decoherence with population decay, and derived the equation of motion for the reduced density matrix 
for the qubit. By numerically solving the equation of motion, we have obtained the time evolution of the 
reduced density matrix. It is shown that, by applying periodic tt pulses, we can suppress the decoherence 
with population decay if the pulse interval is much shorter than the correlation time for qubit-environment 
interactions. 
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